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Abstract 

A novel approach to Riemann-Hilbert problems of particular class is 
introduced. The approach is applicable to problems in which the multi- 
plicative jump is set on a half-line. Such problems are linked to some 
Wiener-Hopf problems motivated by diffraction theory. The new ap- 
proach is based on ordinary differential equations: the Riemann-Hilbert 
problem is reduced to finding a coefficient of an ordinary differential equa- 
tion and solving this equation. The new method leads to an efficient nu- 
merical algorithm and opens a road to new asymptotical and analytical 
advances. 

1 Introduction 

Homogeneous Riemann-Hilbert problems are under consideration. A rigorous 
formulation of these problems and some important results related to them can 
be found in pQ. These problems are constituted in finding a vector function 
analytical (or meromorphic) on the complex plane and having a multiplicative 
jump on a smooth curve. Traditionally, if the curve coincides with the real axis 
the problem is referred to as the Wiener-Hopf factorization problem. Wiener- 
Hopf problems are directly linked with diffraction theory, hydrodynamics and 
some other branches. Numerous important canonical diffraction problems (such 
as diffraction by a half-plane and reflection by waveguide end) can be solved 
by the Wiener-Hopf method [2J. All scalar Wiener-Hopf problems and only 
some particular vector Wiener-Hopf problems can be solved explicitly. The 
solvable problems include a very important Khrapkov's class [3J 0] admitting 
commutative factorization. For other problems some approximate techniques 
exist (see e.g. [5]). All problems studied in the papers mentioned above arise 
from either wave diffraction theory or solid mechanics. 

In the current paper a novel technique is proposed to address the Riemann- 
Hilbert problems possibly not belonging to known factorizable classes. The 
approach is based on the ordinary differential equations. While it docs not 
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provide explicit solutions, it opens some new possibilities for asymptotical and 
analytical studies of the Riemann-Hilbcrt problems. Also it leads to an efficient 
numerical algorithm. 

The current paper continues a research started for the problems of wave 
scattering by a diffraction gratings. A general framework has been developed 
for these problems, comprising application of the so-called embedding formula, 
spectral equation and the problem of reconstruction of an ODE coefficient. The 
current work is dedicated to a possible generalization of the last step of the 
procedure. 



2 The main result 



2.1 Problem formulation 

Consider the following matrix Riemann-Hilbert problem. Let straight contour 
C connect the point b with b + ioo. The complex plane is cut along C. The left 
edge of the cut is denoted by C - , and the right edge is denoted by C + . 

Let there exist an open strip VL of non-zero thickness embracing the cut C, 
such that the point b is located on the boundary of f2 (see Fig. []}. Let 2x2 
matrix function M(z) be analytical in Q, and continuous at z = b. Let also be 



M(z) — > I as z — > b + too, 
where I is a 2 x 2 identity matrix. 
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Figure 1: Contour C and domain Q 



The Riemann-Hilbert problem studied in this paper is formulated as follows. 
It is necessary to find 2x2 matrix function U analytical on the complex plane 
cut along contour C, having no zeros of the determinant on the cut plane, and 
obeying the following condition of the cut: 

U+(z) = ir(z)M(z), (2) 
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where U + (z) and U~ (z) are the limiting values of the function U on the edges of 
the cut. Note that the right multiplication by M is used in problem formulation. 

Additionally, demand that the function U(z) tends to I as \z\ — > oo. Thus, a 
canonical solution is looked for (in Gakhov's sense, [T]). It is assumed also that 
U has no stronger than power growth near z = b. 

2.2 Motivation of the problem 

The problem formulation can be considered as a particular case of a standard 
Riemann-Hilbert problem set on a line going from b + ioo to b — ioo. Namely, 
C is a part of this line, and M can be defined on the rest of the line as a unity 
matrix. However, the author has in mind another link to a standard matrix 
factorization problem. Namely, let b belong to the upper half-plane, and let the 
contour C coincide with the real axis. Let the unknown matrix function U be 
analytical everywhere except C and be discontinuous on the contour C . Let 
the jump on the contour C' be described by the relation 

U(z + Oi) = UO - Oi) A(z), zeC (3) 

where A(z) is a known matrix coefficient. Finally, let A be an algebraic function 
of z defined on the whole complex plane z and having branch points only at 
z = ±b. The last condition is quite strong, however this is not a rare occasion 
in the diffraction-motivated factorization problems. 

Denote by U the continuation of U from the lower half-plane. According to 
(|3"T) this continuation is given by 

U(z) = U(z)A" 1 (z), Im[z] > 0, (4) 

i.e. the function U(z) can be continued onto the complex plane cut along C. 

Denote by the upper indices ± the values of U and A related to the edges 
of the cut C. Obviously, 

U+ = U(z)(A+(z))-\ IT = \J(z)(A-(z))-\ (5) 

Finally, 

fj+(z) =fj-(z)A-(^)(A+(z))- 1 , zeC (6) 

Thus, function U obeys a "half-line" Riemann-Hilbert problem of the type 
introduced above with 

M(*)=A-(z) (A+0))- 1 . 

It is possible to reformulate a lot of practical diffraction problems as half-line 
Riemann-Hilbert problems. 
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2.3 "Ordered exponential" notations 



Let B(t) be a 2 x 2 matrix function analytical in some domain of complex 
plane. Define the ordered exponential of B (denoted by OE[B]) as a solution of 
differential equation 

> j 

_OE[B](r 2l ri) = B(r 2 ) • OE[B](r 2 ,ri) (7) 
with initial condition 

OE[B](n, n ) = I. 

Here t\ and t 2 belong to the domain of analiticity of B(r). Obviously, in this 
case OE[B](t2,ti) does not depend on the contour along which the ordinary 
differential equation is solved. 

If a large area of variation of r 2 is considered, in which B(r 2 ) can have poles, 
it is necessary to indicate a contour 7 connecting t\ and r 2 , along which equation 
is solved. In this case the notation looks like OE[B](7). 

Obvious properties of the ordered exponential are as follows: 

OE[B]( 7l + 72 ) = OE[B]( 72 ) • OE[B]( 7l ), (8) 

OE[B](- 7 ) = (OE[B]( 7 ))- 1 . (9) 

Here and below 7 i + 72 denotes concatenation of the contours, — 7 corresponds 
to the contour 7 passed in the opposite direction. 



2.4 Formulation of the main result of the paper 

Theorem 1 Let r(r) be a function defined and analytical in £1 and continuous 
at t = b such that for all points z 6 (C \ b) the following equation is valid: 

= M(z), (10) 

where contour j z (shown in Fig. [JP lies within Q. z in the left of M0\) plays role 
of a fixed parameter, while t is the independent variable. Then solution of the 
Riemann-Hilbert problem formulated above is given by the formula 

(6,6 + ioo), (11) 

where the ordinary differential equation is solved along contour C. 
Proof 

Obviously, relation (TiT]) defines a function analytical on the complex plane 
cut along contour C. Thus, it is necessary to prove ©. Let be z £ C. Deform 
contour 7z as it is shown in Fig. [3] The contour becomes a concatenation of 7 _ 




U(z) = OE 



r(r) 
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Figure 2: Contours for the OE equation 



going from +ioo to b along the left edge of the cut, and the contour 7 + going 
from b to +ioo along the right edge of the cut. According to ©, ©, and ([TU]) . 



Z — T 

However, according to (|TT|) 
r(r) 



-l 

(-7+) ) OE 



r(r) 



( 7 _) = M(«). 



(12) 



OE 



Z — T 

thus ([2]) is proven □ 



(7-)=U + (z), 



OE 



r(r) 



Z — T 



(~ 7+ )=U-(z), (13) 



3 Interpretation of the main result 

Although the proof of Theorem [T] is simple and straightforward, it would be 
useful to demonstrate another interpretation of this theorem. Let j3 be a point 
on contour C, and let Cp be a part of C going from j3 to -\-ioo. Formulate a 
family of Riemann-Hilbert problems, indexed by /3, having cut Cp instead of C 
(all other conditions of the problem remain the same). The coefficient of each 
problem belonging to the family is a restriction of function M(z) on contour Cp. 
Note that M does not depend on (3. 

Denote by U(/3; z) the family of solutions for these problems. Consider the 
derivative 

WGM = ^U(#z). (14) 

W(/3; z) taken as a function of z is isomonodromical to U(/3; z) in the following 
sense: it is analytical on the plane with the same cut Cp, and the multiplicative 
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Figure 3: Contours 7+ and 7. 



jump on this cut is described by the same coefficient M(z). Moreover, for each 
j3 function W(/3; z) decays as \z\ — > 00, and this function grows no stronger than 
a power function near z — (3. 

The isomonodromy between W and U can be used as follows. One can prove 
that 

W(/3;z) = R(/3;z)U(/3) (15) 

where R(/3; z) is a rational function of z for each fixed j3. To prove this, define 
R as 

R = wr 1 . 

Obviously, R has no singularities on the plane cut along Cp, except maybe the 
point z = j3. The discontinuity of W on Cp is described by the relation 

R+(/3; z) = W+(/3; z) (U+(/3; z))" 1 = W~(fi; z) M(z) M-\z)(U-(fi; z))' 1 = 

W-^^CU-^z))- 1 =BT(p\z), (16) 

i.e. R is continuous on Cp. Thus, R is single valued, has a power singularity at 
z = /?, and decays as \z\ — > 00. Therefore it is rational. A detailed study of the 
singularity at z = /3 shows that R should have form 

W*) = ^- (17) 

Taking into account that the Riemann-Hilbert problem degenerates as /3 — > 
b + ioo, and \j(b + ioo; z) — I, the representation pT|) follows from (|T5|) rewritten 
as a differential equation 

^U(/3;z) = R(/3;z)U(/?;z). (18) 
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Note 1. The method proposed here is based on embedding of the consid- 
ered Riemann-Hilbcrt problem into a family having parameter j3 such that the 
differentiation with respect to the parameter provides an isomonodromy. In this 
context an elementary problem can be considered, for which M is a constant 
matrix with respect to z (or there are several cuts Cj and different constant 
matrices Mj defined on them) . It is not necessary to embed this problem into a 
family of other problems, since differentiation d/dz provides an isomonodromy 
itself. The argument similar to the one provided above shows that U in this 
case is a solution of a Fuchsian equation, and Mj are the monodromy matrices 
of this equation. 

Note 2. Variation of the end point of the cut is not the only possible way 
to organize a family, in which the differentiation with respect to a parameter 
provides an isomonodromy. Another way is to study matrix M(a; z) depending 
on some parameter a and construct a differential operator 

/ n d , . d 
H = ai(a; z)— + a 2 (a; z) — 

such that H[M] = 0. Then it is possible prove that 

H[\J](a;z) = R(a; z) Ufa; z), 
where R is a single- valued function with respect to z. 

4 Examples 

4.1 Scalar example 

Consider a scalar example. Let it be necessary to find function u(z) obeying 
the following restrictions: 

• The function should be regular on the complex plane cut along the line 
C = (b,b + ioo). 

• The values on the edges of the cut should be connected via the relation 

u + (z) =u~(z)m(z), (19) 

where u + are the values on the right edge of the cut, u~ are the values 
on the left edge of the cut, m(z), z e ft is known function tending to 1 at 
infinity. 

• The function should grow no faster than algebraically at b and tend to 1 
at infinity. 

Obtain a "traditional" solution for this problem. Note that a multiplicative 
jump of u corresponds to an additive jump of logu. Thus, the solution can be 
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obtained by using Cauchy's integral: 



logu(z) 



1 

2ni 



b-\-ioo 



logMO) 



(20) 



Now construct a "new" solution for this problem. According to a scalar 
analog of Theorem [TJ it is necessary to find a function r(r) defined in f2, such 
that an analog of (| 10[) is fulfilled: 



OE 



r(r) 



{lz) = m(z), 



This equation can be simplified using the fact that 



OE 



-(r) 



(7 2 ) = exp 



-(r) 



dr 



(21) 



(22) 



log(m(r)) 
2iri 



r(r) 



dr 



(23) 



(24) 



Thus one can find r(r) explicitly: 

r(r) = 

Using the scalar analog of (|TT|) 

m(z) = exp • 

get the solution coinciding with (|20l) . 

4.2 Chebotarev— Khrapkov matrices 

The most important class of matrices admitting explicit factorization is the class 
introduced by Chebotarev |6] and Khrapkov [3]. Later on, these matrices have 
been studied by numerous researchers. Among the most important works there 
are g\EM E2- 

Here the Khrapkov's 2x2 matrices are taken as an example. These matrices 
have the following form: 



M(z) = c(z)I + p(z)L, L(z) 



l{z) m(z) 
n(z) —l{z) 



where l(z), m(z), n(z) are some polynomials such that 

f(z) = l 2 (z) + m{z)n(z) 



(25) 



(26) 



is a polynomial of degree not higher than 2. It is assumed that M — > I as 
z — » b + ioo. 
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It is known that the Khrapkov's matrices can be factorized commutatively, 
but for the price of more flexible behaviour at infinity. Thus, the coefficient 

R(r;,) = ^- 



in (ITT1) and (jTTJJ) should be replaced by a form not necessarily decaying as \z\ — > 
oo: 

R(t;z) = ^I+^L(z), (27) 

Z — T Z — T 

where £ and r\ are some scalar functions. 

The advantage of using the Khrapkov's matrices is their commutativity 
which leads to the closed-form formula for the ordered exponential: 

OE[R(r; z)]( 7 ) = exp{|} (cosh (y/f(z)fj) I + sinh (v?(i)*)) ^L) , (28) 



where 



£ = / -^-dr, f, = / -^-dr, (29) 

J-y Z-T J 1 Z - T 

Thus, in the analog of (JTUJ 

OE[R(r;z)]( 7 ,) = M(z) (30) 

contour j z can be changed to a small closed loop encircling the point z. Ac- 
cording to this, ([50)1 can be solved explicitly: 

{(T) = ilog(c 2 (T)-/(T)p 2 (r)), (31) 



Finally, the solution of the problem is as follows: 

U(z) = ex P {£} ^cosh (V/Rt?) I + sinh (V/(^) -^==^ , (33) 

where 

6+200 6+ioo 
£=- / ^r, ,- = - / J^dr. (34) 



Solution in the form (13"3"|) coincides with the one described in [3]. 
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4.3 A link to Weinstein-class problems 



The current paper can be considered as a generalization of a previous work by 
author dedicated to Weinstein-class problems. Problems of scattering by peri- 
odic diffraction gratings composed of absorbing screens have been studied. The 
first problem having the geometry shown in Fig. @] a) is a classical Weinstein's 
problem [IT] . The period of the system consists of a single screen. The incident 
wave has wavenumber kg (such that k$d ^> 1) and a grazing incidence angle. 
The problem is considered in the parabolic approximation. In |12) this problem 
has been reformulated as the following scalar OE equation: 

OE[r(r)/V^T](7,) - 1 - e lz , (35) 

where 7 Z is the contour introduced above. This equation can be reduced to an 
integral equation with difference kernel and then solved explicitly. It is worth to 
note that this problem has been solved in [IT] by using Wiener-Hopf method. 




a) b) 



Figure 4: Geometry of Weinstein-class problems 

In the second problem of Weinstein's class (see Fig. H]b)) the period of the 
grating contains two screens. In [T3 this problem has been reduced to matrix 
OE equation 

( 7z ) = M(z), (36) 
+ *»). 

The same diffraction problem can be reduced to a Wiener-Hopf problem in 
which the matrix 

K(j3;y) = M(j3-y 2 (d 1 +d 2 )/(2k )) (38) 

should be factorized with respect to variable y for any real value of parameter (3. 

Application of the equation (|36j) enabled one to construct an efficient nu- 
merical solution of the diffraction problem having geometry shown in Fig. |4]b. 
Moreover, using the asymptotic analysis of equation (|3"6")l it became possible to 
prove that the reflection coefficient for this structure tends to — 1 as the angle on 
incidence tends to zero. This result is practically important since it guarantees 
high values of Q-factor in corresponding resonator problems |14j . 



OE 



r(r) 
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5 Solving ( 1101) numerically 



5.1 Parametrization of equation (1101) 

Here equation (jlOp is addressed. It is shown below that the unknown matrix 
r(r) depends on two unknown scalar parameters. These parameters are revealed 
here, and the problem for finding these parameters is formulated in the simplest. 

Let the logarithm of the eigenvalues of M (and the eigenvalues themselves) 
be distinct everywhere on C . Represent M in the form 



M(t) = T(t) 



Ai(t) 
A 2 (r) 



T- X (r) 



(39) 



Here T is the matrix whose columns are eigenvectors of M. Matrix T depends 
on two scalar parameters. Almost everywhere it can be parametrized as follows: 



T(r) 



1 1 

*i(r) Ht) 



(40) 



Values Ai t 2 and ti^C?") are known. 

The left-hand side of (flOl) can be rewritten as follows: 



OE 



r(r) 



( 7z ) = F • OE 



r(r) 



(41) 



where a z is a loop of a small radius e starting and ending at z + ie and encircling 
z in the positive direction, and 



F = OE 



r(r) 



Z — T 

Let r(r) be represented in the form 

' Ci(r) 



(z + ioo, z + ie). 



r(r) = H(r) 









C 2 (r) 



(42) 



As 







OE 



r(r) 



(er 2 ) — > H(z) • exp < — 27ri 



{- 2 "( < f > 6W )}■"■'«»)■ (43) 



It follows from (|4"0)) that the eigenvalues of r(z) are defined by the eigenvalues 
of M(z): 



( h2 (z) = —log(\ h2 (z)). 



(44) 



Thus, to find r one needs only to find H(r). 

Again, the columns of H(t) are the eigenvectors of r. Almost everywhere 
matrix H can be parametrized as follows: 



H(r) 



1 1 

ft x (r) h 2 (r) 



(45) 
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So our aim is to find /ii,2(t). 

For z,f3 £ C , Im[/3] > Im[z] define the function 



V(/3;z) =OE 



r(r) 



(r^) 



where contour r Z)| g is shown in Fig. [5] Obviously, V obeys the equation 

A V (/3; z) = — L- (r(/3) V(/3; z) - V(/3; z) r(/3)) . 



Im[r] 



C 



L z,/3 



Re[r] 



Figure 5: Contour T 



Since V is adjoint to to 



OE 



r(r) 



its eigenvalues are the same, and it can be written as follows: 

V(/3; z) = QG9; z) • ( Al ^ } ^ ) • (Q(/3; z)){z)-\ 

QCM = 



where 



1 1 

qi((3;z) q 2 ((3;z) 



Taking f3 —> z obtain the relation Q(z; z) = H(z), and thus 

91,2(2; z) = h 1:2 (z). 
Taking f3 —> b + ioc obtain Q(oo; z) = M(z), and thus 

91,2(6 + ioo; z) = t 1>2 {z). 



(46) 



(47) 



(48) 
(49) 

(50) 
(51) 
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Elementary calculations demonstrate that equation (jT7|) is equivalent to the 
following system of differential equations: 

Note that these equations are independent from each other. Note also that 
they are Riccati equations. Finally, the reduced problem can be formulated as 
follows: 

For given functions £1,2 (t), £1,2(1") , r G C find functions /ji,2(t), t 6 C 
such that for any z G C equation \52\) possesses a solution qi^(/3,z) on the 
segment f3 G (z,b + ioo) of C obeying boundary conditions H50\). i51\). 

5.2 Numerical algorithm 

The process of solving of the Riemann-Hilbert problem using the new method 
consists of two major stages. On the first stage the coefficient r(r) is found on 
the contour C (it is shown above that it is sufficient to find two scalar functions 
hi^{ T ))- The coefficient is found by solving equation dTU|) or (|52l) . On the 
second stage equation (fTTj) is solved to find function U(z) at some set of points 
Zk- The points at which the U should be computed do not necessarily belong to 
C. Moreover, in diffraction applications C is a half-line parallel to the imaginary 
axis, while it is necessary to find U on the real axis. 

For the first stage contour C is meshed, i.e. a dense array of nodes Tj is 
placed on it. Since the contour is infinite, it is necessary to choose a reasonable 
upper limit B instead of 6 -Moo (see Fig. [5J). The nodes are numbered by indices 
j = 1...N, where 1 corresponds to the node t\ = B (infinity), and index N 
corresponds to "nv = b. 



Im[T,z] 



C 



T. 

J 



points for stage 1 



T N =b 



points for stage 2 



J 



Re[r,z] 



Figure 6: Nodes for numerical computations 



At each tj matrix M(r) is represented in the form (f39|) . (|40[) . i.e. the values 
Ai,2(Tj) and ti^fa) are computed. The values £1,2 (7j) are then computed by 
applying formula (144)) . The aim is to compute the values hi^ijj)- 



13 



For the "infinity" point t\ the following values are assigned: 

fti,2(n) = ti, a (ri). (53) 

This is a natural choice for the asymptotics of the unknown coefficient, since 
M(t) -> I as r -> 6 + ioo. 

Then a loop over j = 2 ... N is performed. At the jth step of the loop 
the values hi^fa) should be computed. Thus, when the jth step is performed 
all values /ji,2(ti) • • • ^l, 2(73-1) are already found. On the jth step of the loop 
equations (j5"2")l are solved on the segment (n, Tj-i) for q± t 2 using Runge-Kutta 4 
method. The initial values are set as 

gi,2(n)=ti,2(7j-) (54) 

following from (|5"Tj) . The values 2(75-1) are found. Then equations (|52|) are 
solved on the segment (73-1, Tj) (only one step is performed) by Euler method. 
This method does not require the values of the right-hand side at the end of the 
segment to be known. The values 91,2 (73) are found. The assignment 

h,2(rj) = 91,2(7-3-) (55) 



is performed, following from (|50|) . 

On the second stage for each Zk the array v(rj)/(zk — Tj) is computed by 
(I52 L (H5 j) . Then the value 



U(z fc ) = OE 



r(r) 



(7-jV,Tl) 



is is found by Runge-Kutta 4 method. 

To demonstrate the feasibility of the new method the Khrapkov's matrix of 
the simplest form has been taken: 

M(z) = I+l(_ 1 z (56) 

The matrix is set on contour C = (2i, ioo). The aim was to find the unknown 
function U at the points Zk belonging to the line lm[z) = 1.8. 

The results obtained by the new method can be compared with the exact 
solution provided by Khrapkov's formulae (|3"3"]) . One has to take into 

account, however, that Khrapkov's formulae lead to the solution Ukh(z) that 
tends to some constant Uoo ^ I as \z\ — > 00. Thus, the Khrapkov's solution 
should be corrected as follows: 

U= (Uoor^UkhXz) (57) 

For the computations the value of B has been chosen to 80«, the step on the 
contour C was equal to 0.02. The results for Re[U] 1 1 and Re[U] 1 2 are shown in 
Fig. El Solid lines represent the Khrapkov's solution corrected by formula ([57)1 . 
dots represent the solution obtained by applying the novel ODE-based method. 
One can see that the agreement between two exact formula and the new method 
is very nice. 
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Figure 7: Results of numerical computations: RefU^ 1 (left), ReflT^ 2 (right). 

6 Conclusion 

The paper presents a new analytical result concerning Riemann-Hilbert problem 
with the multiplicative jump set on a half-line. The solution can be represented 
in the form with the unknown coefficient r(r) obeying equation (1101) . Un- 
fortunately, the new result does not lead to explicit solving of a general matrix 
factorization problem, but it provides an efficient numerical technique and gives 
some hopes to achieve further progress in asymptotical and analytical studies. 
The efficiency and robustness of the numerical method is linked with the fact 
that the equation (|10[) has the "Volterra" structure, i.e. the unknown coefficient 
can be found point- by-point. The asymptotical advantages of the new method 
has already been demonstrated for similar problems of Weinstein type. Finally, 
some analytical progress, hopefully, can be achieved by studying the analyti- 
cal structure of the auxiliary function r(r) for the case of rational or algebraic 
matrix M. 

The new method brings the Riemann-Hilbert problem back into the context 
of the ordinary differential equations (from which it initially emerged) , but in 
a slightly unusual form. Instead of studying the derivative with respect to the 
main variable z the derivative with respect to the parameter /3 is considered. 

A new method is proposed for Riemann-Hilbert problems belonging to a 
particular class. As it is shown in the paper, the class naturally arises as a 
reformulation of a problem of factorization of a matrix A on the real axis if 
the matrix is given as an algebraic function having one branch point in the 
upper half-plane (or in the lower half-plane). Such case is quite common in the 
diffraction theory, but it restricts the class of problems to those having only one 
propagation speed. If there exist several wave modes with different speeds the 
new methods needs to be generalized to admit several half-infinite contours Cj . 
Such a generalization can be obtained by allowing the coefficient R to have 
several pole singularities. 

The paper considers only matrices 2x2, however most of the results can be 
easily converted for the case of matrix of arbitrary dimension. 
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